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Abstract In situ measurements of [OH], [HO 2 ] (square brackets denote species concentra- 
tions), and other chemical species were made in the tropical upper troposphere (TUT). [OH] 
showed a robust correlation with solar zenith angle. Beyond this dependence, however, [OH] 
did not correlate to its primary source, the product of [0 3 ] and [H 2 0] ([0 3 ]*[H 2 0]), or its sink 
[NO y ], This suggests that [OH] is heavily buffered in the TUT. One important 
exception to this result is found in regions with very low [0 3 ], [NO], and [NO y ]. 
Under these conditions, [OH] is highly suppressed, pointing to the critical role of NO 
in sustaining OH in the TUT and the possibility of low [OH] over the western Pacific 
warm pool due to strong marine convections bringing NO-poor air to the TUT. In 
contrast to [OH], [HO x ] ([OH] + [H0 2 ]) correlated reasonably well with [0 3 ]*[H 2 0]/[N0 y ], 
suggesting that [0 3 ]-[H 2 0] and [NO y ] are the significant source and sink, respectively, of 
[HOx]- 
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1 Introduction 

The hydroxyl radical (OH) is the most important oxidant of a number of short-lived 
stratospheric ozone-depleting substances in the tropical upper troposphere (TUT) 
(Brioude et al. 2010). Understanding its distribution in the TUT is thus critical for 
predicting the ozone depletion potential of these very short-lived substances. The OH 
distribution in the TUT is also important for controlling the delivery of sulfur to the 
stratosphere and, thereby, the amount of stratospheric aerosol (Weisenstein et al. 1997; 
Rae et al. 2007). There are few OH measurements currently available in the tropical 
upper troposphere, particularly in regions of active troposphere-to-stratosphere 
exchange. 

In the TUT, OH is produced from 0 3 and H 2 0 photochemical reactions (Kley et al. 1996; 
Wennberg et al. 1998; Jaegle et al. 2001): 


O 3 H- h'v — » O 2 + 0( D) 

(Rl) 

0( 1 D) + H 2 0 — ► 20H 

(R2) 


Due to uplift associated with strong convection, the composition of air in the TUT, 
particularly in the tropical Pacific region, is greatly affected by near-surface air (Kley et al. 
1996; Folkins et al. 2006). As a consequence, 0 3 concentrations in the TUT can reach very low 
levels (<20 ppb) in the Western Pacific warm pool region. Furthermore, TUT H 2 0 concen- 
trations over the warm pool are very low due to the very low ambient temperatures (Newell 
and Gould-Stewart 1981). Since O3 and H 2 0 are source species for OH, low concentrations of 
these species may lead to low [OH] (square brackets denote species concentrations 
hereafter) and some model studies appear to support this conjecture (e.g. Kley et al. 
1996). Other model studies (Crawford et al. 1999) showed a very weak dependence of 
OH on [0 3 ] and [H 2 0]. This weak dependence was due to a strong HO x ([OH] + 
[H0 2 ]) source from autocatalytic production through methane oxidation (Wennberg 
et al. 1998; Jaegle et al. 2001). 

Observationally, [OH] is found to be largely independent of [0 3 ] and [H 2 0] in both the 
lower stratosphere (LS) (Wennberg et al. 1994b; Hanisco et al. 2001) and free troposphere 
(Rohrer and Berresheim 2006). In both cases, [OH] values were found to be proportional to 
ambient photolysis rates; namely, a linear combination of photolysis rates for O('D) 
and CH 2 0 in the LS and the photolysis rate for 0(*D) alone in the free troposphere. 
These proportionalities suggest that either the sinks for OH are correlated with OH 
sources, or that OH is highly buffered in these regions. We note that the latter is a 
sufficient (but not necessary) condition for the near constant global OH levels 
reported by Montzka et al. (2011). 

In this study, we show OH measurements from an extended set of high-altitude aircraft 
flights that provide a unique dataset for examining the role of OH in the TUT. We show that 
[OH] in the TUT behaves very much the same as [OH] in the lower stratosphere as reported by 
Hanisco et al. (2001). While the product of [O 3 ] and [H 2 0] varies by about 3 orders 
of magnitude, most [OH] data deviate less than a factor of 2 from a functional fit (a 
linear combination of photolysis rates for O(’D) and CH 2 0) of stratospheric [OH] vs. 
solar zenith angle (SZA) (Hanisco et al. 2001). The largest anomaly from this near- 
constancy is the association of very low [OH] values with air masses depleted in 
[NO], 
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2 Measurements 

The dataset used for this work was collected using a suite of in situ instruments onboard the 
NASA ER-2 high-altitude research aircraft during the NASA Stratospheric Tracers of 
Atmospheric Transport (STRAT) field mission (1995-1996). OH and HO 2 were 
measured with a laser- induced fluorescence instrument (Wennberg et al. 1994a). O3 
was measured by a dual-beam UV photometer (Proffitt and McLaughlin 1983). A 
Lyman-a hygrometer was used for H 2 0 measurements (Weinstock et al. 1994). 
Measurements of NO and NO y were made with a chemiluminescence detector (Gao et al. 
1997). CO was measured using a laser absorption spectrometer (Webster et al. 1994). Temper- 
ature and pressure were obtained using the onboard meteorological measurement system (Scott 
et al. 1990). 

Latitudes sampled by ER-2 flight tracks during STRAT ranged from 2°S to 60°N. 
Plights during STRAT originated either from Moffett Pield, CA or Barbers Point 
Naval Air Station (now Kalaeloa Airport), HI. Only flights based at Barbers Point 
were used in this work (1 to 23°N). A significant amount of flight time (~76 %) was 
spent in the lower stratosphere during these flights. Only upper tropospheric data 
between 9 km and the tropopause with 0 3 mixing ratio below 150 parts per billion 
(ppb) are used for this study. Tropopause heights along the flight tracks were 
determined using temperature profiles taken by the Microwave Temperature Profiler 
instrument (Denning et al. 1989). 

Both NO and NO y are tightly correlated with 0 3 , while H 2 0 is clearly dependent on 
ambient temperature (Pig. 1). The N0 y -0 3 correlation has a slope of 0.005, approximately 
twice that observed in the tropical LS in this dataset and 1.6 times that seen in the higher- 
latitude LS (Fahey et al. 1996). The CO mixing ratio spans a relatively small range (70± 
20 ppb) and is loosely anti-correlated with 0 3 . As Fig. la and b suggest, NO and NO y in this 
dataset are also highly correlated (slope = 0.41, r 2 =0.82). [N0 2 ] values during daytime (SZA 
< ~85°) were in a photochemical steady state with [NO] via reactions R3 and R4 (Gao et al. 
1 997), and were much lower than [NO] values because formation from [NO] via its reaction 
with ozone is suppressed due to the low [0 3 ] and temperature values while photolysis of N0 2 
remains efficient: 


NO -b O 3 — > N0 2 A 0 2 


(R3) 


N0 2 + hw -> NO + O (R4) 

[NO x ] ([NO] + [N0 2 ]) is essentially the same as [NO] shown in Fig. lb. Air parcels with 
very low [0 3 ] and [NO y ] comparable or close to those reported by Kley et al. (1996) (Fig. la 
inset) were encountered. As in Kley et al. (1996), these air parcels likely experienced recent 
convection and wet removal of water-soluble species prior to sampling. The absence of 
enhanced NO levels implies that this convection was not accompanied by large amounts of 
lightning (Skamarock et al. 2003). 


3 Results and discussion 

[OH] and [HO x ] are shown as functions of SZA in Fig. 2. [H0 2 ] is much greater than [OH] 
and, therefore, is the dominant component of [HO x ] (Wennberg et al. 1998). In Fig. 2a and b, a 
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Fig. 1 NO, NO y , H 2 0 and CO 
mixing ratios shown as functions 
of O3. The dataset includes 
approximately 3,000 points, each 
of which is a 10-s average 
corresponding to approximately 
2 km of flight track. Latitude, 
longitude, temperature, and alti- 
tude are shown as color scales in 
panels a, b, c, and d. Linear fits of 
NO and NO y are shown in panels a 
and b as black lines 
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Fig. 2 [OH] and [HO x ] shown as 
functions of solar zenith angle 
(SZA). The data points correspond 
to those in Fig. 1. A functional fit 
to OH concentrations measured in 
the lower stratosphere (Hanisco 
et al. 2001) is shown in panels a 
and b as the solid lines. The dotted 
and dashed curves are the Hanisco 
fit ±50 % as visual aids 
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functional fit to [OH] in the lower stratosphere (Hanisco et al. 2001) is shown as the black 
curve. This fit is expressed as: 

[OH(SZA)] = ci Jj + c 2 J 2 

J = J 0 exp[-£,(sec(x)-l)] m 

Jj : Jo = 8 x 10~ 5 , t = 2.2, x = 0.85 x SZA v ’ 

h : Jo = 6.4 x 1(T 5 , C = 0.67, x = 0.86 x SZA 

where Jo is in s -1 , C! = 1.5x to 10 molecules cm -3 s, and C2=1.75x 10 10 molecules cnf 3 s. 

A few measurements made by the ER-2 while sampling its own exhaust plume (Hanisco 
et al. 1997) are also shown in Fig. 2 as large black circles (to be discussed below). 

The TUT [OH] values in Fig. 2a and b essentially follow the same relationship with respect 
to SZA as those measured in the lower stratosphere (Hanisco et al. 2001). The primary 
production rate of OH via reactions R1 and R2 is proportional to the product of 0 3 and 
H 2 0 concentrations ([0 3 ]*[H 2 0]). In this dataset, where [0 3 ]-[H 2 0] varies over 3 orders of 
magnitude, [OH] does not show any correlation to [0 3 ] , [H 2 0] or [NO y ], The absence of 
correlation of [OH] to its primary production source (R1 and R2) suggests that low [0 3 ] and 
[H 2 0] do not necessarily lead to low [OH]; consequently the primary production rates alone 
are insufficient to predict [OH], 

The lack of correlation of [OH] to its source gases is not entirely surprising. [OH] and 
[H0 2 ] are tightly coupled via reactions R5-R8: 


OH + CO(+0 2 ) -> H0 2 + C0 2 

(R5) 

OH + 0 3 — > H0 2 + 0 2 

(R6) 

H0 2 + NO — » OH + N0 2 

(R7) 

H0 2 + 0 3 — > OH + 20 2 

(R8) 


The cycling time between OH and HO 2 is only a few seconds, while the lifetime of HO x is 
much longer (on the order of 30 min) (Jaegle et al. 2001). [OH] and [H0 2 ], therefore, are in a 
photochemical steady state under virtually any sunlit condition. [OH] is thus controlled not 
only by its production reactions, but by the production reactions for HO 2 , the loss reactions of 
[HO x ] and the redistribution between [OH] and [H0 2 ] via reactions R5 to R8 (McKeen et al. 
1997; Crawford et al. 1999; Jaegle et al. 2001). [HO x ] is controlled only by its sources and 
sinks. The combination ofRl andR2 is one of the mainHO x sources, andNO y is a major HO x 
sink (Wennberg et al. 1998; Jaegle et al. 2001; Hanisco et al. 2001). In this dataset, [0 3 ]*[H 2 0] 
showed no correlation with [NO y ] (r 2 =0.17). This lack of correlation is expected because 
[H 2 0], [0 3 ], and [NO y ] are independently determined by an air parcel’s temperature history, 
photochemical history, and the combination of photochemical history and wet removal 
processes, respectively. [HO x ] showed little or no correlation with [0 3 ]*[H 2 0] (r 2 =0.12; data 
with SZA >65° are excluded due to their apparent correlation with the changing solar 
radiation or SZA) and l/[NO y ] (r 2 =0.49). This is not surprising since any possible single 
correlation between HO x and one of its sources or sinks can easily be obscured by variations in 
other sources and sinks. [HO x ] did show a significant correlation with [0 3 ]*[H 2 0]/[N0 y ] 
(Fig. 3, / =0.72 ), consistent with the notion of [0 3 ]*[H 2 0] (via R1 and R2) and [NO y ] (via R9 
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Fig. 3 Correlation plot of HO x density (10-s average) and [0 3 ]*[H 2 0]/[N0 y ]. Data points with solar zenith 
angles greater than 65° are excluded to minimize the effect of changing solar radiation 


and RIO) being a significant source and sink, respectively. This correlation can be seen in 
Fig. 2c and d, where at any SZA less than 65°, low [HO x ] corresponds to low [ 03 ]*[H 2 0 ] and 
high [NO y ], and vice versa. Because neither [0 3 ]*[H 2 0] nor [NO y ] are the sole source and sink 
for [F10 x ], a tight correlation between [HO x ] and [0 3 ]-[H 2 0]/[N0 y ] is not expected. 

The theoretical [H0 2 ]-to-[OH] ratio can be calculated using reactions R5-R8 and measured 
[0 3 ], [CO], [NO], pressure, and temperature: 

[HO 2 ] _ k OH+C o[C O] + [0 3 ] 

[OF1] £//o 2 +vo[NO] + £//0 2 +o 3 [O 3 ] 

where £ x denotes the respective reaction rate. The calculated [H0 2 ]-to-[OH] ratio is compared 
with measurements, expressed as a ratio of the calculated ratio to the measured ratio, in Fig. 4. 
As shown, the calculations generally agree with the measurements to within 30 %. The NO 
signal-to-noise ratio is low at extremely low and high SZAs (Fig. 2b), and hence results at SZA 
< 8.5° and >78° are not included in the plot. The contribution of R 6 is negligible and R7 
dominates the F10 2 reactions. The steady-state approximation used in deriving Eq. (2) is 
reasonably accurate even at the lowest [NO] values, in which case the F10 2 self-reaction (R12) 
is still approximately a factor 3 slower than reactions R7 and R 8 . 

Reactions involving HNO 3 and HNO 4 , significant components of NO y , are sinks of F10 x 
The key reactions are: 

ITNO 3 + OFI — > FI 2 0 + N0 3 (R9) 


HN0 4 + 0H^H 2 0 + N0 2 + 0 2 (RIO) 

R7, R9, and RIO are a set of buffering reactions for OFI when R9 and RIO are major sinks 
for F10 x and when NO is correlated with F 1 N 03 and HNO 4 : Fligher [HNO 3 ] and [FINO 4 ] lead 
to lower [HO x ], but the associated higher [NO] forces a lower [F10 2 ]-to-[OFI] ratio in the F10 x 
partitioning. An example where this buffering effect breaks down can be seen in the ER-2 
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Fig. 4 Ratio of 10-s averages of measured and calculated H02-to-OH ratios are shown as a function of solar 
zenith angle (SZA) 


plume sampling events (large black circles in Fig. 2); here NO produced by the ER-2 engine 
(Fahey et al. 1995) clearly forced [OFI] out of its nominal relationship with SZA while [F10 x ] 
values show no apparent response to this additional [NO]. (Note that OH production may be 
enhanced by HONO in the ER-2 exhaust (Hanisco et al. 1997), but the photolysis of HONO is 
too slow to affect the [H0 2 ]-to-[OH] ratio.) When NO y is too low to be the main HO x sink, this 
buffering effect will also break down. 

When [NO y ] is low, [HO x ] is limited by the quadratic self-reactions: 


OH + H0 2 -> H 2 0 + 0 2 

(RU) 

HO 2 + HO 2 — ► H 2 O 2 + O 2 

(R12) 

OH + H 2 0 2 -> H 2 0 + H0 2 

(R13) 


and reactions with CH 3 OOH (Jaegle et al. 2001). These self-reactions partially explain the 
relatively small changes in [HO x ] when [0 3 ]*[H 2 0] varies over 3 orders of magnitude. 
Another possible reason for the small HO x variations is additional, more constant HO x sources 
via ketone photochemistry (Me Keen et al. 1997) and autocatalytic production through meth- 
ane oxidation (Crawford et al. 1999; Jaegle et al. 2001). 

While [0 3 ] and [H 2 0] show little effect on [OH], low [NO] and [NO y ] conditions may lead 
to anomalously low [OH] due to the breakdown of the N0-HN0 3 -HN0 4 buffer. Linder these 
conditions, while NO y is too low to be the major sink of HO x , the HO x values still are limited 
due to its self-reactions. While HO x values are self-limited, the low values of NO may not be 
able to recycle sufficient OH from H0 2 (R7). As shown in Fig. 2, very low [OH] values (more 
than a factor of 2 lower than the functional fit of Hanisco et al. (2001)) were observed during 
the STRAT mission (5°-10° SZA in Fig. 2a and b). These low [OH] values were associated 
with low [0 3 ] (<2 x 10 1 1 enf 3 ), [NO] (<1 x 10 8 cm -3 ), and [NO y ] (<2 x 10 8 cm -3 ). Note that the 
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corresponding [HO x ] (i.e. HO 2 ) was actually at its highest within the STRAT dataset. In this 
case the recycling rate of HO 2 back to OH was too slow to keep up with the loss of OH 
through reaction with H0 2 , resulting in the decreased [OH], This link between low [NO] and 
low [OH], if characteristic of air lifted to the TUT over the western Pacific warn pool via 
convection, will have important implications. The western Pacific warm pool is known to be a 
region of large-scale deep convection, which brings a large amount of 03 -poor air to the TUT. 
It appears that [NO] is also likely to be low in these air parcels both from limited STRAT data 
and the fact that lightning is rare in marine convection (Crawford et al. 1997; Christian et al. 
2003; Cecil et al. 2012). These low values occur in a region where air is transported into 
stratosphere (Schoeberl et al. 2013). Low [OH] values may affect the breakdown rates of 
various ozone depleting species (Brioude et al. 2010). We note that [OH] may be enhanced in 
TUT air parcels originating from land-based convection, where lightning is common (Cecil 
et al. 2012 ). 


4 Conclusions 

OH, H0 2 , and other chemical species were measured in the TUT region during the STRAT 
aircraft campaign. While [OH] values show a robust relationship with SZA, they are essen- 
tially independent of [0 3 ] and [H 2 0]. The majority of measured [OH] values, therefore, can be 
expressed using the simple formula (Eq. 1) reported by Hanisco et al. (2001) for stratospheric 
measurements. 

The robust [OH]-SZA relationship is likely due to strong buffering. Awell-known buffering 
mechanism (NO favoring OH within the partitioning of HO x , and HNO 3 and HNO 4 limiting 
HO x concentration) appears to play a role when NO y is high. However, it is unclear that 
whether this mechanism alone is sufficient to account for the robust [OH]-SZA relationship. 

Low-[OH] air parcels that fall well outside the [OH]-SZA relationship were encountered at 
times during Pacific tropical flights. These air masses are characterized by very low [O 3 ], 
[NO], and [NO y ], We postulate that the low [NO y ] conditions render the NO-HNO 3 -HNO 4 
buffer ineffective, while low [NO] slows down the recycling of OH from H0 2 , thus leading to 
low [OH] in these air parcels. This potential link between low [NO] and low [OH] has 
important consequences such as lowering the removal efficiency of surface pollutants over 
the west Pacific warm pool and allowing more efficient transport of very short-lived ozone- 
depleting substances and other compounds such as S0 2 into the stratosphere. 
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